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ABSTRACT 

We present a population synthesis study of the observed properties of the magnetars 
investigating the hypothesis that they are drawn from a population of progenitors 
that are more massive than those of the normal radio pulsars. We assume that the 
anomalous X-ray emission is caused by the decay of a toroidal or tangled up field that 
does not partake in the spin down of the star. Our model assumes that the magnetic 
flux of the neutron star is distributed as a Gaussian in the logarithm about a mean 



value that is described by a power law <& = $o 



9Mr 



G cm2 (8M0 s^ Mp ^ 45M0) 



where Mp is the mass of the progenitor. We find that we can explain the observed 
properties of the magnetars for a model with $0 = 2 x 10^^ G cm^ and 7 = 5 if 
we suitably parametrise the time evolution of the anomalous X-ray luminosity as an 
exponentially decaying function of time. Our modelling suggests that magnetars arise 
from stars in the high mass end (2OM0 ^ Mp ^ 45 Mq) of this distribution. The lower 
mass progenitors are assumed to give rise to the radio pulsars. 

The high value of 7 can be interpreted in one of two ways. It may indicate that 
the magnetic flux distribution on the main sequence is a strong function of mass 
and that this is reflected in the magnetic fluxes of the neutron stars that form from 
this mass range (the fossil field hypothesis). The recent evidence for magnetic fluxes 
similar to those of the magnetars in a high fraction (~ 25%) of massive 0-type stars 
lends support to such a hypothesis. Another possibility is that the spin of the neutron 
star is a strong function of the progenitor mass, and it is only for stars that are more 
massive than ~ 2OM0 that magnetar-type fields can be generated by the a—cu dynamo 
mechanism (the convective dynamo hypothesis). In either interpretation, it has to be 
assumed that all or a subset of stars in the mass range ^ 2OM0 — 45M0, which on 
standard stellar evolution models lead to black holes via the formation of a fall-back 
disc, must give rise to magnetars. 

Unlike with the radio pulsars, the magnetars only weakly constrain the birth spin 
period, due to their rapid spin-down. Our model predicts a birthrate of ^ 1.5 — 3 x 
10~^ yr~^ for the magnetars. 

Key words: pulsars: general, stars: neutron, stars: early-type, stars: magnetic fields 



1 INTRODUCTION 

The anomalous X-ray pulsars (AXPs) and Soft Gamma Re- 
peaters (SGRs) have spin down properties indicative of mag- 
netic fields of ~ 10^''"^® G which places them at the high end 
of the neutron star field distribution. The source of the qui- 
escent X-ray luminosity in both the SGRs and the AXPs is 
unclear, but it is generally believed that it is associated with 
the decay of some component of the magnetic field. 

The standard model of iDuncan fc ThompsonI (119921 ) 
assumes that the fields are generated at the time of for- 
mation of the neutron star by an efficient a — uj dynamo 
that operates at low Rossby numbers and requires millisec- 
ond birth periods. Neutron stars born with initial periods 



Pi are predicted to generate large scale magnetic fields of 
3 X 10^^ G (1 ms/Pi) under optimum conditions, which is 
more than adequate to explain the fields in magnetars. 

A consequence of the rapid initial spin predicted for 
these models is that the supernova explosions that create 
the magnetars are expected to be an order of magnitude 
more energetic than ordinary core-collapse supernovae, if 
one makes the standard assumption that angular momen- 
tum is lost by magnetic braking and not by gravitational 
radiation or due to emission in a jet. However, the en- 
ergetics of some well studied supernova remnants associ- 
ated with magnetars, appear to suggest that their formation 
may not always be accompanied by a hypernova explosion 
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IjVink fc Kuipenl200d ) . The origin of magnetic fields in mag- 
netars therefore remains unclear. 

Recent observations have provided two additional clues 
on the origin of magnetars. Firstly, several magnetars have 
been associated with massive star clusters. The magne- 
tar candidate CXO J164710. 2-455216 h as been associated 
with the open cluster Westerlund 1 (JMuno et al.l |2006| ) 
and S GR 1806-20 wit h the giant Galactic HII complex 
W31 IjFiger et al.ll2005h . both suggesting progenitors with 
M~4OM0. Massive (~20M(T) ) progenitors appe ar also to 
be indica ted for SGR 1627-41 ("C orbel et ahll 19991 ) and SGR 
0526-66 (jKlose et al.l |2004). Interestingly, the progenitors 
have masses in the range ~ 20 — 45M0 where interac- 
tion with a fall-back disc is expected to play an important 
role in determining the ultimate outcome of stellar evolu- 
tion - nam ely whether a neutron star or a black hole forms 
JHeger et al. 2003). 

A second clue comes from the discovery that more mas- 
sive stars tend to have higher magnetic fiuxes than less mas- 
sive stars, and that the incidence of magnetism on the main 
sequence increases with mass. Some 25% of O and early B 
stars in the Orion cluster are strongly magn et ic with fields 
in the range 600 - 1500 G (|Petit et all l|2008t ). lDonati et all 
1120021 )') so that the progenitors of the magnetars are also 
likely to be strongly magnetic. Fossil magnetic fluxes similar 
to those observed in the magnetars may already be present 
in stellar cores prior to collapse and could give rise to global 
magnetic fields of ~ 10^^ G either directly through magnetic 
flux conservation, or indirectly by acting as seed fields for 
a dynamo. In this context, we note that while the dynamo 
mechanism is effective in generating weak flelds of the order 
of the equi-partition valu e in the convecti ve cores of stars 
with radiative envelopes (JBrun et al.ll2005l ). much stronger 
flelds can be generated by the dynamo mechanism if a fossil 
field is present in the radiative region (iF eatherstone et al.l 
l2006l ). 

In this paper, we carry out a population synthesis study 
to analyse the observed numbers of magnetars and their 
magnetic field and period distributions allowing for X-ray 
selection effects. Our basic premise is that the birth mag- 
netic field of a neutron star is related primarily to the mass 
of the progenitor star Mp. Our aim is to determine the form 
of this relationship in the hope that it may provide some 
insights into the origin of magnetars. 



2 THE POPULATION SYNTHESIS MODEL 

2.1 Modelling the spin down of magnetars 

We assume that the magnetic fi eld of a magnetar is com- 
posed of two components, as in IPons et al.l (|2007|): an un- 
derlying dipolar component Bd and a crustal component 
Be- We assume that the dipolar component is a global field 
that does not decay on a time scale of ~ 10* yr, the lower 
limit to the field decay time scales in normal radio pulsars 
(e.g. see the most recen t population sy ntheses calculations of 
IVranesevic et al.l 1120041) .iFerrario fc Wickramasing he (2006) 
(hereafter FW), iFaucher-Giguere fc Kaspil (|2006l )) and is 
solely responsible for the spin down of the magnetar which 
occurs along constant field lines. The crustal component is 
envisioned to be a toroidal or tangled field which does not 



partake in the spin down of the star, decays on time scales 
of ^ 10'* — 10^ yr, and is the source of the anomalous X-ray 
and 7— ray emission (see section [3] for a discussion of the al- 
ternative hypothesis that the spin down is due to a decaying 
field). 

The apparent upper limit of ~ 12 s for the spin 
periods of magnetars has to be attributed to the ex- 
pected rapid decline in the X-ray luminosity predicted by 
field decay models (e.g. see th e period clustering study of 
iGolpi. Geppert. fc Page! l|2000f )). Models of neutron stars 
with a heat source in the crust predict that cooling will 
occur initially at a nearly constant luminosity after which 
the luminosity will dec ay exponentially to negligible values 
(jKaminker et al.l 120061 ). but there are no detailed models 
that investigate the dependence of decay rates on initial 
magnetic field and envelope composition. In the absence of 
such calculations, we have adopted the following parametri- 
sation to describe the decline of the crustal luminosity 



where 



L = Lo 



Td ~ TdO 



Bd 



10" G 



10" G 



Bd 



t 

Td 



yr- 



Here, Lq and Tdo are constants, and we have assumed that 
the birth crustal luminosity is proportional to the intrinsic 
dipolar field. We have also allowed for the possibility that 
the decay time is a function of the magnetar's field strength 
Bd via the parameter 5. 



2.2 Parametrisation of magnetic field 

Recently, there have been detailed stellar evolution calcula- 
tions for stars in the mass range I2M0 ^ Mp ^ 35M0 that 
have attempted to incorporate magnetic fields and rotation 
(iHeger et al.ll2005l ) . These calculations allowed for fields gen- 
erated in the differentially rotating radiative regions of stars 
and their effect on the transport of angular momentum, on 
the assumption that a dynamo mechanism is in operation. 
However, the basic premise on which these calcu lations were 
made has recently been seriously challenged bv lZahn et al.l 
(2007 ). Therefore, the role played by magnetic fields and 
rotation in the evolution of stars still remains uncertain. 

It is likely that the magnetic flux of the progenitor star, 
^p, and to a lesser extent the initial angular momentum, 
both play an important role in determining the magnetic 
flux <I>NS of the nascent neutron star. In our modelling, we 
assume that the progenitor mass Mp is the primary variable 
that describes the birth magnetic flux of the neutron star. 

Our modelling procedure is as discussed in FW, except 
for our treatment of the birth magnetic flux and birth spin 
period. Briefly, we begin with a dynamical model for the 
Galaxy, an initial mass function, and an initial-flnal mass 
relationship for neutron stars. We then follow the spin evolu- 
tion and motion of neutron stars and calculate the properties 
of the currently observed population of magnetars allowing 
for X-ray selection effects. 

In the present study, we assume that the birth magnetic 
flux of the neutron star is distributed as a Gaussian about 
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a mean that is given by a power law 



$NS = *o 



/ Mp 



V9Mg 



G cm 



(8M0 s: Mp s; 45M(: 



(1) 

with dispersion criogSpjg. For consistency with our previous 
study of radio pulsars, we assume log $0 = 25.2 and deduce 
the exponent 7 from the study of magnetars. The birth spin 
period is assumed to be given by a Gaussian distribution 
with mean Po and dispersion apg. 

The birth magnetic field of the neutron stars is calcu- 
lated from 



Bns 



•I- 



7r(i?is 



(2) 



where 7?ns is the radius of the neutron star in cm given by 
the meiss-radius relationship for neutron stars which depends 
on the equation of state (see FW for further details). 



2.3 Allowing for X-ray selection in the ROSAT 
survey 

We have chosen to use the subset of magnetars discovered 
from the ROSAT All-Sky survey to constrain our models. 
In order to achieve this, we need to account for interstellar 
extinction and X-ray selection effects. 

We have implemented the methodology laid out by 
iGill fc He vl (2007) for finding the percentage of the magne- 
tars in our calculations that yield a ROSAT limiting count 
rate of 0.015 counts s~^ in the 0.1 — 2.4 keV energy range 
JHiinsch et al.lll999l l. 

Our synthetic objects have been assigned a two- 
component model (power law and bla ckbody) modulate d 
by interstellar absorption. Hence, as in lGill fc Hevll ((20071), 
the number of photons in the energy range E — E + dE is 
given by 



N{E) dE = 



m 



f ^ 



VlkeV 



+ 



V2 



8.1E^ 



fe4T4(e-SA2-_i) 



-NHa{E)^j^ 



(3) 



Here, r/i is the power-law normalisation in units of pho- 
tons keV~^cm~^s~^, a is the power law energy index, kT is 
the temperature in keV and 772 is the blackbody normalisa- 
tion in units of L^g/Dlo, where Z/39 is the luminosity in units 
of 10^^ erg s~^ and Dio the distance in units of 10 kpc. The 
neutral hydrogen co l umn density Nh is obtained using the 
[Foster fc Routledgel (120031 ) model of the galactic HI distri- 
butio n. T he photoelectric absorption cross-sections cr{E) are 
from iBal ucihska-Church & McCammon (1992) used with 
the ISM abundances of iWilms et al.l (|2000l ) . In order to con- 
vert the absolute photon fiuxes into the photon counts of the 
X-Ray telescope (XRT) detector, we have used the ROSAT 
calibration guidelines and the XRT effective area combined 
with the PSPC sensitivity. 

In our modelling, the magnetars have been assigned a 
temperature and an energy index according to Gaussian dis- 
tributions with means and spreads kT = 0.5, ut ~ 0.1 and 
a = 3.5, (Jq — 0.5 respectively. The relative contribution of 



the black body and power law components to the total X- 
ray luminosity has been randomly assigned wit h values con- 
sistent with the observed distribution given bv lPerna et al.l 
(2001 ). 



3 RESULTS AND DISCUSSION 

We have computed a series of models for different values M* , 
7, (Tiogso, ^0, o-po, Lq, Tdo and S until a model is obtained 
that satisfies the following criteria: (a) the model must enve- 
lope and provide a reasonable representation of the observed 
period and magnetic field distributions and (b) the total 
number of observed magnetars must be less than the theo- 
retically predicted number after X-ray selection effects are 
taken into consideration. Given the small number of magne- 
tars and the unknown number of active, but currently dor- 
mant magnetars (e.g. objects such as SGR 1627-41 that has 
only recentl y re-activated after a nearly 10 year long quies- 
cent period IjPalmer et al.ll2008l )). the real Galactic distribu- 
tions is not well enough sampled to justify a more detailed 
modelling. 

The total number of active magnetars in the Galaxy 
depends strongly on M* because of the rapidly declining 
initial mass function. Our calculations show that in addition 
to the initial mass function, the parameters 7, $0, o"iog#o, 
Lo, Tdo and 5 together determine the magnetic field, period, 
and period derivative distributions. 

Our results are rather insensitive to Po and apg be- 
cause of the rapid initial spin down of magnetars. Hence, 
this study cannot exclude the possibility that magnetars 
are bor n at millisecond periods as r equired by some recent 
models IjGeppert fc Reinhardtll2006l ). or that the birth spin 
period may be a function of magnetic field as assumed in 
FW. For simplicity, in our models we have adopted Pq = 3 
seconds and apg = 1 second, although the FW prescrip- 
tion with an upper initial period cutoff at 8 seconds is also 
acceptable. 

We present in Figure 1 our favoured model with M* — 
2OM0, log$o = 25.2, criog*o = 0.4, Lo = 10^^ r^o = 
5 X 10^ yr, and S = 1.3. The corresponding birthrate is 
3 X 10~^ yr~^. This model best envelopes the field, period, 
period derivative, characteristic age, luminosity and tem- 
perature distributions and predicts 26 active magnetars af- 
ter applying the ROSAT detectability criteria. We exclude 
the possibility of a significantly lower Ad* (e.g. ^ 18Mq) 
on the grounds that it will lead to a conflict with the birth 
rate of radio pulsars which we assume arise from stars with 
A'l ^ M* . On the other hand, the use of a significantly 
higher A4* (> 22Mq) leads to too few magnetars and incon- 
sistencies with the observed distributions. 

We have also constructed a P — P diagram which 
is shown in Figure 2 where our ROSAT selected objects 
are compared with the observations of all magnetars with 
known P. Thus, among the observed objects we have also 
included SGR 0526-66, which is located in the LMC, and 
CXOU J010043- 721134 which is in the SMC, since their po- 
sition in the P — P diagram is unlikely to depend on which 
galaxy they reside. 

The predicted number of active magnetars from our 
favoured model can be reconciled with the observed number 
of 5 ROSAT detected sources if we assume that only ~ 20% 
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Figure 1. Our model (solid line) overlapped to the magnetars data sample (dashed histograms) taken from 
http://www. physics. mcgill.ca~pulsar/magnctar/main. html and with the addition of the recently confirmed member of the AXP 
class, 1E1547.0-5408 (=PSR J1550-5418) llCamilo et al.ll2007l'l . 



of magnetars were in an active state which brought them 
to a high enough luminosity to be detected by this satel- 
lite. In this case, the actual birthrate of magnetars would be 
~ 3 X 10~^yr~^, which of course depends on our basic as- 
sumption that all stars in the mass range 20 — 45Mq produce 
magnetars. However, we note that the total number of ob- 
served magnetars in the Galaxy is only 14. If no additional 
magnetars exist in our Galaxy, that is if all magnetars have 
already been discovered, this means that we have produced 
too many magnetars from progenitors in the mass range 
20 — 45M0 . This may suggest that some 50% of stars in this 
mass range will produce black holes. If this is the case, the 
birth rate of magnetars would reduce to ~ 1.5 x 10"'^ yr"'^- 

Studies of SGR bursts conducted bv iKouveliotou et al.l 

([199J) indicated that t here cannot be more tha n 7 active 
SGRs in the Galaxy and IKouveliotou et al.l (|199a ) suggested 
that magnetars are born at a rate of about 0.1 per cen- 
tury. The recent study carried out bv lGill fc Hevll (|2007l ) was 
based on a study of the five AXPs detected in the ROSAT 
All-Sky Survey, and indicated a birthrate of 0.22 per century 
with their progenitors being massive main sequence stars. 
Both these estimates are generally consistent with the birth 
rate that w e dedu ce. Another recent study carried out by 
iMuno et al.l l|2008l ) who searched for magnetars in archival 
Chandra and XMM-Newton observations of the Galactic 
plane, yielded a birthrate in the range 3 x 10~^ — 6 x 10~^ 
yr~^. The upper end of this range is excluded by our calcu- 
lations. 

On the fossil field hypothesis, the high value indicated 
for the magnetic flux index 7 may simply reflect the intrinsic 
magnetic flux distribution on the main sequence. 

Observations of main sequence stars show that the max- 



imum magnetic fluxes observed on the main sequence map 
on to the magnetic fluxes of the highest fleld magnetic white 
dwarfs (10^-10^ G ) and neutron stars (10"-10^^ G) rather 
well. However, due to sensitivity limitations of polarimetric 
observations, only the upper end of the magnetic field distri- 
bution of main sequence stars (the strongly magnetic stars) 
can be observed, so we only have a partial picture of mag- 
netism on the main sequence. Unfortunately, most stars in 
the mass range 8 — 2OM0 that give rise to radio pulsars have 
magnetic fiuxes well below the currently observed range on 
the Main Sequence so that the index 7 cannot be empirically 
estimated. We note however that magnetic fluxes similar to 
those inferred in magnetars occur in ~ 25% of massive B 
and 0-type stars. 

If we somehow dismiss the close similarities between the 
magnetic fluxes of massive main sequence stars and magne- 
tars by putting them down to mere coincidence, alternatives 
to the fossil fleld hypothesis need to be explored. It is possi- 
ble that the mass of the progenitor determines the spin of the 
nascent neutron star and thereby the strength of a dynamo 
generated field. Support for this h ypothesis comes from the 
calculations of iHeger et al.l (|2005l ) which allow for angular 
momentum transport by magnetic fields generated by dif- 
ferential rotation during stellar evolution. These show that 
more massive stars tend to produce more rapidly spinning 
neutron stars. For a 35M0 progenitor (the highest mass they 
considered), they predict that the neutron star will have a 
spin period of 3 ms which is rapid enough for the generation 
of a a magnetar-type field by the a — uo dynamo mecha- 
nism of iDuncan fc ThompsonI l| 19921 ). However, for this to 
be a viable explanation of the magnetars, the a ~ lu dy- 
namo mechanism would need to be effective for neutron stars 
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log(P) 

Figure 2. Filled squares: observed magnetars 

(http: //www. physics. mcgill.car^pulsar /magnetar/main. ht ml) 
with the addition of 1E1547.0-5408 fCamilo et all ' 20 07|); filled 
circles: magnetars as derived from our modelling (see text). The 
solid line is an empirically determined boundary (a "magnetar 
death line") given by log(P) = 8.41og(P) - 20. 



that arise from progenitors with significantly lower masses 
(~ 20 — 22M0) or the derived birthrate of magnetars would 
be too low to explain the observations. Furthermore, since 
the calculationsof Hcgcr ct al. (2005) have been seriously 
challenged bv lZahn et al.l (|2007l ). it appears that we are still 
waiting for more stellar evolution calculations, that allow for 
both magnetic field and rotation, to show us whether or not 
millisecond rotation periods can result from massive stars. 
If such periods can be achieved, this would give support to 
the idea that the dynamo model is a viable alternative to 
the fossil field hypothesis, although it is likely that a fossil 
field could still play the role of a seed field. 

It is relevant that we comment on the alternative to 
our basic hypothesis on spin evolution, namely that the field 
that decays is also the field that drives spin evolution dur- 
ing the magnetar ph ase. We have considered this p ossibility 
using "Avenue C" of lColpi. Geppert. fc Pagj (|2000l ) , but ex- 
cluded it from our present considerations because it failed 
to populate the high field end of the observed magnetar field 
distribution. It is conceivable that with the use of different 
field decay parameters for avenue C we may also be able to 
model the data, but we expect that such a model will need 
to have field characteristics during the magnetar phase that 
are so similar to those adopted in the present calculations 
that our major conclusions will remain largely unchanged. 

In conclusion, we note that the origin of the crust al field 
compo nent in magnetars remains unresolved. IZahn et al.l 
l|2007l ) have shown that the shearing of a poloidal field of 
fossil origin in a differentially rotating radiative region can 
lead to the generation of a toroidal field so that field com- 
plexity of the type that appears to be required in magnetars 
may arise naturally from stellar evolution, when the com- 
plex interplay between rotation and fossil magnetic fields is 
taken into consideration. On the other hand, since fall-back 
discs are expected to play a role in the evolution of stars 
in the mass range '^ 20 — 45M0, it is tempting to specu- 



late that the crustal field component that characterises the 
magnetars may have its origin in the interaction of the fall- 
back disc with the magnetic field of the stellar core, in cases 
where a magnetar is the outcome of stellar evolution. If this 
is the case, one would expect that most of the mass in the 
fall-back disc would be magnetically ejected. 
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